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We present calculations of the magnetic ground state of CS2CUCI4 in an applied magnetic field, with the aim 
of understanding the commensurately ordered state that has been discovered in recent experiments. This layered 
material is a realization of a Heisenberg antiferromagnet on an anisotropic triangular lattice. Its behavior in a 
magnetic field depends on field orientation, because of weak Dzyaloshinskii-Moriya interactions. We study 
the system by mapping the spin-1/2 Heisenberg Hamiltonian onto a Bose gas with hard core repulsion. This 
Bose gas is dilute, and calculations are controlled, close to the saturation field. We find a zero-temperature 
transition between incommensurate and commensurate phases as longitudinal field strength is varied, but only 
incommensurate order in a transverse field. Results for both field orientations are consistent with experiment. 

PACS numbers: 75.10.Jm, 75.25,+z, 75.40.Cx 
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I. INTRODUCTION 

Experiments on the spin 5=1/2 triangular lattice antifer- 
romagnet CS2C11CI4 have mapped out its properties in great 
detail over the last few yearsii^^i^ It has attracted particular 
interest because its low-dimensionality, frustrated interactions 
and small spin are all features expected to promote quantum 
fluctuations. It is therefore a system for which long-standing 
proposals 7 of quantum disordered states may be relevant at 
intermediate energy scales, despite the fact that it has conven- 
tional Neel order at the lowest temperatures. Experimental 
studies are facilitated by the small strength of exchange inter- 
actions, which makes it possible to reach the saturation mag- 
netization with accessible fields and also results in relatively 
low excitation energies. Unusual features of the excitations 
revealed in inelastic neutron scattering 4 have attracted consid- 
erable theoretical interest. 8 i 9 i 10 i n i 12 i 13 i 14 

The magnetic phase diagram of CS2CUCI4 as a function of 
temperature, magnetic field strength and field orientation is re- 
markably rich: it reflects the competition between exchange, 
Dzyaloshinskii-Moriya (DM) and Zeeman energies. In zero 
field, incommensurate spiral long range order is observed at 
temperatures below Tm — 0.62K. The ordered moments lie 
in an easy plane (the b-c plane, see Fig. selected by the 
DM interactions. An applied magnetic field has different ef- 
fects depending on whether it is transverse to this plane (along 
the a-axis) or longitudinal (within the b-c plane). In a trans- 
verse field at low temperature, ordered moments cant out of 
the easy plane, lying on a cone, and an incommensurately or- 
dered phase is observed for all field strengths up to the satu- 
ration field of B^ r = 8.51T. The magnetization, incommen- 
surate ordering wavevector and transverse order parameter in 
this phase have, we believe, been well accounted for in earlier 
theoretical work by the present authors in collaboration with 
Coldea, using the 1/5 expansion, together with dilute Bose 
gas techniques close to the saturation field. 15 Thermal effects 
in the dilute Bose gas have been examined recently in Refs. 5 
andd 

By contrast, behavior in longitudinal fields is less well un- 
derstood. From a theoretical viewpoint, competition between 



magnetic field and easy plane anisotropy adds complexity 
to the problem. Using a classical (large 5) calculation, we 
found in earlier work 15 two incommensurate low-temperature 
phases: a distorted cycloid spin structure in low fields, and a 
cone state at higher fields, with a ferromagnetic ally aligned 
state above the saturation field. Initial neutron diffraction 
studies in a magnetic field directed along the c-axis identified 
a distorted cycloid for fields below B c < 1.4T and elliptical 
incommensurate order for 1.4T < B c < 2. IT. 17 Subsequent 
measurements detected a further incommensurate phase in a 
narrow range just below the saturation field B c cr = 8.02T, 
for 7. IT < B < B c cr (Refll5l). and most recently, interven- 
ing commensurate order, for 2. IT < B c < 7. IT (Ref0). 
All incommensurate phases have magnetic Bragg peaks at the 
field-dependent wave vectors Q(B) — (0, Q(B), 0), whereas 
in the commensurate phases Bragg peaks are observed on the 
Brillouin zone boundary, at Mi, M2 = (0, |, ±5). 

In the present paper, we examine theoretically the zero tem- 
perature phase diagram of CS2CUCI4 in a magnetic field with 
a strength close to the saturation value. We focus on the ob- 
served differences in the effects of longitudinal and transverse 
fields. Behavior close to the saturation field is of special inter- 
est because in this regime quantum fluctuations are under the- 
oretical control, even for 5 = 1/2. Reversed spins form a gas 
of hard core bosons which is dilute, and interactions between 
bosons can be treated exactly at low density by summing two- 
body scattering processes to obtain an effective interparticle 
potential»A&I2i22i2L Using this approach we have shown previ- 
ously 15 that the system has an incommensurate cone ground 
state at fields just below the saturation value, for both field 
orientations. In the following we add to this with the finding 
that the dilute Bose gas treatment yields a first order quantum 
phase transition to a commensurately ordered ground state be- 
low a longitudinal field strength of B^ l — 0.96B^ r . Since 
commensurate order appears at a field strength Bt\ which is 
close to the saturation field B c cr , it takes place under condi- 
tions for which our calculations are expected to be reliable. 
For a system in a transverse field the same method gives only 
an incommensurate cone state. Results for both field orienta- 
tions are therefore consistent with experiment. Our ideas can 
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be further tested experimentally by comparing the magnetic 
structure we predict for the commensurate phase with obser- 
vations. 



II. MODEL AND CALCULATIONS 

To proceed, we first recall the Hamiltonian appropriate for 
CS2CUCI4 , which has been determined with precision by in- 
elastic neutron scattering measurements of the energies of sin- 
gle spin-flip excitations from the fully polarized, high field 
state. 3 Magnetic moments in a single layer of the material lie 
on sites of an anisotropic triangular lattice as shown in Fig.^ 
The dominant exchange interaction J = 4.34K is along the 
crystallographic b axis. A weaker exchange J' = 0.34J 
acts on the zig-zag bonds, and a much weaker interaction 
J" = 0.045 J couples layers. DM interactions are symmetry- 
allowed on the zig-zag bonds. The spin Hamiltonian can be 
written as H = Ho +Hdm + Kb, where Ho, Hdm and Kb 
include the isotropic, DM and Zeeman interactions, respec- 
tively. Here 

Ho = E JsSr ■ Sr +< 5, (1) 

R,<5 

where S denotes bond vectors connecting neighboring sites 
and Jg represents J, J' or J", as shown in Fig. ^ The DM 
interaction is 

Hdm = ' Sr. x [Sr+5! . (2) 

R 

The vector D± = (±D, 0, 0) (with D = 0.053 J) is perpen- 
dicular to the layers and alternates in sign between even and 
odd layers because they are inverted versions of one another. 
The Zeeman interaction is Hb — — X)r H • Sr, where H is 
a reduced Zeeman field with components H a = g a ^,sB a . 

This spin model can be treated as a lattice gas of hard core 
bosonsii£ii2i£2i2ii22 In the dilute limit, all quantum effects are 
incorporated exactly by summing ladder diagrams for the in- 
teraction verticesi^ 3 ^ A Hartree Fock treatment of these ef- 
fective interactions can be used to obtain the ground state spin 
structure. We describe further details of the calculation sep- 
arately for the two cases of longitudinal and transverse mag- 
netic fields. 



A. Longitudinal Fields 

For definiteness, we consider a field in the c direction. We 
introduce boson creation and annihilation operators ^> R and 
R and set S R = \ - <0 R , S+ = 5 R + zS R = R and 
5 R = 5 R — iS^ = </> R . For the spin commutation relations 
to be satisfied, a hard core constraint on boson number must 
be imposed by including an on-site repulsion U in the Hamil- 
tonian and taking the limit U — > 00. 

The DM interaction is represented by terms cubic in bo- 
son creation and annihilation operators, because for this field 
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FIG. 1 : (Color Online) (A) Magnetic sites and exchange couplings of 
CS2C11CI4 : Even (color) and odd (gray) layers are stacked along the 
a-direction. The signs C5 , refer to the direction of the DM vector 
between zig-zag bonds. (B) Reciprocal space plane (k, I), with the 
projected Brillouin zone, projections of the Brillouin zone center V, 
the incommensurate (Q, — Q) and the commensurate (L, Mi, M2) 
ordering wavevectors indicated. (C) Predicted antiferromagnetic 
state (Eq. [8} with Bragg peaks at Mi and M2, showing the ori- 
entation of ordered moments projected onto the a-b spin plane. Note 
that here, of necessity, spin coordinates are rotated relative to lattice 
coordinates. 



orientation it couples transverse and longitudinal spin com- 
ponents. Close to the saturation field we are concerned ul- 
timately with an effective description involving only bosons 
near minima of the dispersion relation generated by the ex- 
change coupling. As argued previously, 15 momentum conser- 
vation precludes cubic terms in an effective low energy Hamil- 
tonian. Instead, within second order perturbation theory the 
DM interactions generate quadratic and quartic couplings be- 
tween low energy bosons, with magnitude 0(D 2 /J). Since 
\D/J\ <C 1, we neglect altogether the effect of DM interac- 
tions in longitudinal fields. 

At this point a simplification is possible. In the absence 
of DM interactions the distinction between odd and even lay- 
ers disappears. The size of the unit cell in the a-direction is 
halved and that of the Brillouin zone is doubled. We take 
Fourier transforms using this reduced unit cell. The Hamilto- 
nian (omitting a constant) is 

k k,k'.q 

(3) 

where = Jk — J m in is the boson kinetic energy, Jk = 
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1 /2 J2s Jse is the Fourier transform of the exchange cou- 
plings, /i = H c cr — H c is the boson chemical potential, 
and H£ r = Jo — J m in is the saturation field. The num- 
ber of lattice sites is N and the interaction vertex is given 
by V q = 2J q + 2U. Writing k = (h,k,l) as short- 
hand for k = 2n(h/a,k/b ) l/c), the dispersion relation ek 
has two minima located at the incommensurate wavevectors 
k = ±Q = ±(1, l/2 + e ,0), where sin 7re = J'/(2J). On 
reverting to the standard description with two layers in a unit 
cell, these wavevectors are denoted by Q = ±(0, 1/2 + eo, 0) 
and amplitudes are staggered on the two layers. 3 Note that the 
degeneracy of the two minima follows from symmetry of the 
full Hamiltoniar. 15 and is not lifted by DM interactions. 

At this point we can apply standard techniques developed 
for the interacting Bose gas. 23 With negative fj,, the ground 
state is the boson vacuum. Equivalently, with H c > H£ r 
the spin system is fully polarized. With positive /i, the bo- 
son density is non-zero. Condensation of these bosons at low 
temperature is equivalent to magnetic order, and we introduce 
the complex order parameter ipq = (4>q) / >/N. In the dilute 
regime, the full scattering vertex can be obtained by summing 
pair interactions in the particle-particle channel. We denote 
this vertex, for the scattering amplitude of incoming bosons 
with momenta k and k' to outgoing states with momenta k— q 
and k' + q, by T q k,k'- At low density it satisfies the Bethe- 
Salpeter equation 



r 



q,k,k'=Vq - 



V, 



q-q' 



£ k+q' 



£k - £k' 



q',k,k' 



(4) 

This integral equation can be reduced to a set of linear equa- 
tions which are readily solved numerically (see Ref.fl^). 

We now proceed to an analysis of the energy of various 
candidate ground states. This energy can be written in the 
form of a Landau expansion, in terms of one or more order 
parameters tpq. We first review our earlier results, for behavior 
immediately below the saturation field. 15 In this regime one 
expects condensation only at one or both minima of the boson 
dispersion relation, with wavevectors k = ±Q. Allowing for 
two possible order parameters, the energy per site is 



(5) 



with Ti = r 0: Q,Q = r ,-Q,-Q and T 2 = IV-q.q + 
r2Q,-Q,Q- (Here cq = 0, but we retain it for future refer- 
ence). Minimizing the energy with respect to the order param- 
eters, a single component state is favored for Ti < IV It is the 
cone state with spin structure (Sr.) = (y/nQ cos(±Q ■ R + 
a), % /n Q sin(±Q-R+a), 5-riq), where n Q = (/i-e Q )/r 1 
is the condensate density and a is an arbitrary phase. A two- 
component state with equal amplitudes for both order param- 
eters is favored for Ti > T 2 - In this state, ordered mo- 
ments for a spin fan, with (Sr) = (2y / ng cos(Q • R + 
a) cos(/3), 2V" q cos(Q -R + a) sm(/3), 1/2 -4n Q cos 2 (Q • 
R + a)), where tiq = n_ Q = {p — e Cl )/(Y 1 + T 2 ) and a, (3 
are arbitrary phases. A numerical evaluation of Eq. |4]yields 



i?/iV = (e Q - M )(|^ Q | 2 + |^-Q| : 

+ ^r 1 (|^ Q | 4 + |^ Q | 4 )+r 2 ^ Q | 2 |V-Q| 2 



Ti = 3.79 J and T 2 = 4.67 J, so the cone state is selected. It 
has energy per site E/N = — (eq — ^i) 2 /2ri. 

At this stage, we have established that the ground state for 
H = + is the cone state. Next we examine whether there is 
a transition to a commensurate ground state at a larger value 
of /i. A first possible commensurate state is one with order at 
the wavevector L = (1, 1/2, 0), using our description with a 
single layer per unit cell (equivalent to (0, 1/2, 0) in the stan- 
dard notation with two layers per unit cell). This is a natural 
choice because L lies close to Q, and because the wavevector 
of incommensurate order is known to move towards L with in- 
creasing We find nevertheless, repeating the calculations 
we have outlined, but with L in place of Q and with el > 0, 
that the commensurate spin cone and spin fan states are higher 
in energy than the incommensurate states, as shown in Fig. [2] 

An second possible commensurate state is one with con- 
densates at the wavevectors Mi = (l, i, |) and M2 = 
(l, —5, |), in our notation (or (0, |, 5) and (0, — |, |) in 
the standard notation), and it is in fact at those positions that 
magnetic Bragg peaks are observed experimentally. 6 The dif- 
ference AM = M2 — Mi is half a reciprocal lattice vector 
so that umklapp scattering between the condensates is allowed 
and may reduce the energy of the state. We obtain an energy 
per site 



E/N=(e Ml -n) (Wj' + I^/W 

+ \ T A (l^Mj | 4 + |^M 2 | 4 ) + T B \lpM 2 | 2 I^Mi I 



+ (^MiV^ +C-C 



(6) 



where 



0,Mi,Mi 



ro,M 2 ,M 2 — 3.36J, T B 



ro,Mi,M 2 + ^AM.Mi ,Mj — and T c — T^m,Mi,Mi — 

1.30 J. The amplitude Tb is zero because the denominator 
in Eq. (0}, £k+q' + £k'-q' — £k — e k r , depends on only two 
components of q' for k = Mi and k' = M2. This leads to 
vanishing interactions, as is familiar from the example of the 
two dimensional Bose gasi^ 5 . 

The energy of this commensurate state is minimized by set- 
ting t/> Ml = ±#m 2 , giving 



N ~ 



T A 



(7) 



B 



We find that it is lower than the energy of the incommensu- 
rate cone state for fields smaller than = 0.958H^ r : see 
Fig- The proximity of to H£ r justifies our use of the 
low density approximation. The components of the ordered 
moment in this commensurate state are 



R/ 



(Si 



R/ 



n c (cos(Mi • R+ a) ± sin(M 2 • R + a)) , 
(sin(Mi • R + a) =F cos(M 2 • R + a)) , 

- 2n n , 



(8) 



where n c = (fi — e^^/i^A + — Tc) ar, d a is an arbi- 
trary phase. This structure is illustrated in Fig.^ In this com- 
mensurate state neighboring spins are antiferromagnetically 
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aligned along the chains and in adjacent layers, but perpendic- 
ular in adjacent chains. The ± sign in Eq.|S]reflects the two 
distinct ways to arrange this perpendicular orientation. The 
experimental phase diagram in longitudinal field is summa- 
rized in Fig. |2j together with the theoretical results obtained 
here for behavior close to the saturation field, and earlier ones 
for large S. 
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FIG. 2: (Color Online) Ground states in longitudinal fields. Top: 
Energies of states considered [Eqns. J5) and JfJ] as a function of lon- 
gitudinal magnetic field B c , computed using a c = 2.30^ The thin, 
thick, short-dashed and long-dashed lines represent, respectively: the 
incommensurate cone state; the commensurate state with conden- 
sates at Mi and M2; the incommensurate spin-fan; and the commen- 
surate state with condensate at k = L. The incommensurate cone 
state is the ground state for B^x < B c < B c cr , with B c cr = 7.95T 
and B^i = 7.615T. The commensurate state at Mi and M2 is the 
ground state for B c < B^ ri . Bottom: Phase diagram in longitudinal 
fields: (A) from experiment; (B) from dilute Bose gas calculation; 
(C) from large S-expansior.' 5 . Here S denotes an incommensurate 
phase, and AF and AF* indicate commensurate phases, observed 
experimentally to have magnetic Bragg peaks at Mi and M2. 



where Slk = Jk — Jo + H a , evaluated with J" = 0, Z?k = 
—D (sin(k ■ <5i) + sin(k • <S 2 )) and J£ = J" cos(k • S 3 ). Di- 
agonalization of this Hamiltonian yields a dispersion rela- 
tion with two branches eir. Because of DM interactions, the 



branches lack inversion symmetry: ej^ 7^ e^k- 5 Below the sat- 
uration field, order appears at the incommensurate wavevec- 
tors that minimize and e^*££ It is likely that any further 
transition to a commensurate phase will be disfavored by the 
reduced symmetry of the dispersion relation in a transverse 
field. To test this, we have carried out calculations of the ener- 
gies of commensurate states in a transverse field analogous to 
the ones described for a longitudinal field. To simplify these 
calculations, we have treated a model in which the DM vector 
is not staggered between layers, so that the unit cell contains 
only a single layer, rather than two. (See Ref. 5 for an alter- 
native approach.) A similar calculation to the one performed 
in longitudinal fields reveals that the incommensurate state is 
lowest in energy in transverse fields (see Fig. [3ji in agreement 
with experimental data and the 1 / S expansion. 
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FIG. 3: (Color Online) Ground state in transverse fields. Energies of 
various states as a function of transverse magnetic field B a (normal- 
ized with respect to the critical field). The thin, thick, short-dashed 
and long-dashed lines represent, respectively: the incommensurate 
cone state; the commensurate state with condensates at Mi and M2; 
the incommensurate spin-fan; and the commensurate state with con- 
densate at k = L. 



B. Transverse fields 



III. SUMMARY 



In a transverse field, the DM vector is parallel to the field 
direction and so the DM interaction is quadratic in boson 
creation and annihilation operators. Using the standard unit 
cell containing two layers, and introducing separate species 

of bosons for the even and odd layers, the quadratic Hamilto- 

;„3 





~n k -j"+D k j\ 




0eU 




ilk — J" — Dk 







In summary, we have studied the magnetic ground states 
of CS2CUCI4 in the vicinity of the saturation field by treat- 
ing reversed spins as a dilute Bose gas. In a transverse field 
we find only an incommensurate ground state , in agreement 
with experiment and calculations based on the 1/5 expansion. 
By contrast, in a longitudinal field there is a transition be- 
tween an incommensurate state close to the saturation field 
and a commensurate state at lower field. We propose that this 
commensurate state is the one recently observed in neutron 
diffraction^ 
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